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Recent Results from RHIC & Some Lessons for Cosmic-Ray Physicists 

Spencer R. Klein^ * 

'^Nuclear Science Division, Lawrence Berkeley National Laboratory, Berkeley, CA, 94720, USA 

The Relativistic Heavy Ion Collider (RHIC) studies nuclear matter under a variety of conditions. Cold nuclear 
matter is probed with deuteron-gold collisions, while hot nuclear matter (possibly a quark-gluon plasma (QGP)) 
is created in heavy-ion collisions. The distribution of spin in polarized nucleons is measured with polarized proton 
collisions, and photoproduction is studied using the photons that accompany heavy nuclei. 

The deuteron-gold data shows less forward particle production than would be expected from a superposition 
of pp collisions, as expected due to saturation/shadowing. Particle production in A A collisions is well described 
by a model of an expanding fireball in thermal equilibrium. Strong hydrodynamic flow and jet quenching shows 
that the the produced matter interacts very strongly. These phenomena are consistent with new non-perturbative 
interactions near the transition temperature to the QGP. 

This writeup will discuss these results, and their implications for cosmic-ray physicists. 
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Since the first collisions in 2001, RHIC has pro- 
duced a wealth of data on pp, dAn, AuAu and 
CuCu collisions at center of mass energies from 
20 to 200 GeV per nucleon pair. RHIC was built 
to study aspects of QCD; the main foci have been 
the study of cold nuclear matter via dA colli- 
sion, and of hot nuclear matter (the Quark Gluon 
Plasma (QGP)?), via heavy-ion collisions. RHIC 
also studies polarized proton collisions to measure 
the polarized parton distributions, and photonu- 
clear interactions. This data is also of interest for 
modelling cosmic-ray air showers. 

RHIC hosts 5 experiments. There are two large 
experiments, PHENIX and STAR, and 3 smaller 
ones: PHOBOS, BRAHMS and pp2pp. 

PHENIX has two large central spectrometers, 
and a forward muon system [T]. It is optimized 
for particle identification, particularly for leptons. 

STAR is optimized for global event studies, 
with a large acceptance time projection chamber 
for charged particles and a calorimeter to detect 
neutral particles [2] . 

PHOBOS detects charged particles over a very 
large pseudorapidity range, \r]\ < 5.4, and has two 
small spectrometer arms for tracking Pseu- 
dorapidity 77 = — In [tan(0/2)], where 9 is the an- 
gle between the particle direction and the beam 
axis. BRAHMS consists of central and forward 
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spectrometers, with precise particle tracking and 
identificaton in a small solid angle [4^ . pp2pp con- 
sists of Roman pots which track protons scattered 
at small angles; it studies pp diffraction [5]. 

This writeup will review the different RHIC 
physics, starting with polarized proton collisions 
and photoproduction, before moving on to discuss 
cold and hot nuclear matter. 

1. pp and Polarized Proton Collisions 

pp collisions serve two functions at RHIC: tests 
of pQCD calculations with unpolarized pp colli- 
sions, and measurements of the polarized parton 
distributions. 

Jet and single particle cross section data at 
RHIC are in good agreement with recent pertur- 
bative QCD calculations [6|7j . Figure [T] compares 
pQCD calculations with data for 7r° cross sec- 
tions; the agreement is good. This good agree- 
ment is possible thanks to recent next-to leading 
order pQCD calculations and improved fragmen- 
tation functions. 

Polarized proton collisions are used to study 
the spin structure of the nucleon [8]. The total 
nucleon spin is the sum of the quark spins, gluon 
spins and orbital angular momentum within the 
proton. Experiments at SLAG and CERN in the 
19 70s- 1990s showed that the quarks carry a rel- 
atively small fraction of the nucleon spin; this is 
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Figure 1. PHENIX data on 7r° production in 
in 200 GeV pp collisions, compared with NLO 
pQCD calculations plus two different fragmenta- 
tion functions [B]. 



sometimes called the 'spin crisis.' RHIC will mea- 
sure the gluon polarization, via quark-gluon and 
gluon-gluon interactions. The polarized gluon 
distributions are found by comparing the normal- 
ized rates for collisions where the proton spins are 
pointing in the same vs. in opposite directions, 
Although much more data is needed, already it is 
possible to exclude the most extreme models of 
gluon polarization [7]. 

2. Photoproduction 

Photoproduction has been studied in gold-gold 
and deuteron-gold collisions. Virtual photons 
from the electromagnetic field of one nucleus in- 
teract with the other nucleus. Photoproduction 
of the J /ip is sensitive to the gluon distributions 
of the target nucleus [9]. p and 7r"'"7r~7r"'"7r~ final 
states have also been studied [TU] . 

Vector meson photoproduction can occur two 
ways: nucleus 1 can emit a photon which inter- 
acts with nucleus 2, or vice- versa. The two chan- 
nels are indistinguishable, so they interfere. Go- 
ing from one emitter to the other is a parity in- 
version, and vector mesons are negative parity, so 



the two amplitude subtract. At mid-rapidity 

a(6)=ao(6)(l-cos(pT-fo)). (1) 

Here, CTq is the cross section without interfer- 
ence the pt is that of the vector meson, and h 
is the impact parameter (distance between the 
two ion centers at closest approach), h is not an 
observable, so the total cross section is the in- 
tegral of Eq. 1. This interference has been ob- 
served through a reduction in cross section for 

PT < n/{b) m- 

Because of the strong nuclear fields, some re- 
actions can uniquely be studied in heavy-ion col- 
lisons. At RHIC, multi-photon interactions in- 
volving a single ion pair i.e. Au + Au —^ Au* + 
Au* + p via 3 photon exchange (one for each nu- 
clear exchange, and a 3rd to produce the p") has 
also been studied ^Oj. The LHC will reach ^p 
center of mass energies up to several hundred 
GeV, well beyond the reach of HERA [12]. 

3. Cold Nuclear Matter 

By comparing pp and dA collisions the effects 
of the nuclear environment may be studied. The 
nuclear environment is expected to alter parton 
densities, possibly leading to qualitatively new 
behavior [T3]. Gluons have a virtuality [Q^) de- 
pendent transverse area ^{fic)^ / Q"^ . At high den- 
sities (i.e. at low x), partons will recombine, in 
reactions like g+g g. This recombination mod- 
erates the growth in gluon density as x decreases. 
Because of the higher parton density, recombina- 
tion is more significant in nuclei, and comparisons 
of parton distributions in protons and in nuclei 
are sensitive to these changes; the difference mea- 
sured in lepton-nucleon interactions is known as 
shadowing. 

Figure [2] shows the effect of this density on 
particle production. It compares charged-particle 
production in dAu and pp collisions at different 
pseudorapidity [Hj. Higher pseudorapidity cor- 
responds to higher- a; partons from the deuteron, 
and lower-x partons from the gold; as expected 
from recombination, shadowing is larger at lower- 

X. 

Different theoretical studies have used BFKL 
and/or DGLAP evolution, or taken advantage of 
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Figure 2. The ratio of charged particle produc- 
tion in dA and pp colhsions (normahzed to the 
number of nucleon-nucleon coUisions) at different 
pseudorapidity. In the absence of nuclear effects, 
RdAu should equal 1. jTi|. 



the relationship between shadowing and diffrac- 
tion. One interesting new approach, the colored 
glass condensate (CGC) treats the gluon fields in 
the nucleus as a classical field. 

The CGC makes some interesting predictions 
|15j . a CGC should interact coherently as a sin- 
gle object, producing some new effects. When 
a parton interacts with a CGC, the CGC recoils 
coherently; because of it's high mass, the recoil 
is muted, leading to apparent 'monojets'. With- 
out a CGC, 2 parton 2 parton reactions like 
g + g g + g produce two azimuthally back-to- 
back jets. However, if the gluon strikes a CGC, 
the heavy object will recoil slowly, and the recoil- 
ing gluon will produce the single visible jet. This 
process is studied experimentally via 2-particle 
correlations. Figure [3] shows the azimuthal angle 
correlations between a tt*^ produced in the for- 
ward direction and a charged hadron produced 
near mid-rapidity, for pp and dA collisions. Di- 
jet events should produce back-to-back correla- 
tions. These correlations are smaller for dA col- 
lisions than for pp collisions, as expected from 
CGC models. The suppression rises at smaller 7r° 
energies, as expected from a CGC. 

4. Hot Nuclear Matter 

A main experimental goal of RHIC is to search 
for the quark gluon plasma (QGP), an interact- 
ing system of partons in equilibrium [17j . Lattice 
gauge theory predicts that a QGP is produced 
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Figure 3. Azimuthal correlations between a for- 
ward tt'^ and a more central charged hadron in 
(left) pp and (right) dA collisions for (top) 25 < 
< 30 GeV and (bottom) 30 < < 55 GeV 
[16] . The lower (blue) area shows the background, 
while the higher (red) area follows a fit to the data 
(points). The 'S' values give the size of the peak, 
smaller in dA collisions. The correlations are ex- 
pected in dijet events. 



when nuclear matter at low baryochemical poten- 
tial (baryon density) is heated to a temperature 
above about 170 MeV [18]. 

RHIC heats nuclear matter by colliding heavy 
ions. The collisions occur in several stages. Ini- 
tially, the ions collide and their partons inter- 
act, producing new partons. These partons them- 
selves interact, producing still more partons. As 
these interactions occur, the system expands and 
cools. Eventually, the partons form hadrons 
which may themselves interact. Interactions con- 
tinue until the system cools enough that no fur- 
ther inelastic processes are possible; this transi- 
tion is known as chemical freezeout. Later, the 
hadrons separate enough that no further elastic 
collisions occur; this is thermal freezeout. 

Much RHIC data is analyzed in terms of cen- 
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Figure 4. dN/drj for 3 — 6% copper-copper col- 
lisions and 35 — 40% gold-gold collisions. The 
two datasets have very similar Npart, and almost 
identical dN/drj T9]. 



trality, or impact parameter (b). Several meth- 
ods are used to characterize centrality: Npart is 
the number of nucleons participating in the col- 



lision; for a head-on collision, N, 



part 



2A. N, 



bin 



is the number of binary (nucleon-nucleon) colli- 
sions; this is relevant for comparing pQCD par- 
ticiple production between pp and AA collisions. 
For a given b, Npart and N{,in are determined us- 
ing a Glauber calculation. Centrality is also given 
in percentages, such as the — 10% most central 
collisions. 



4.1. Overall Event Structure 

We begin by considering the particle produc- 
tion. Figure |4] compares the pseudorapidity (77) 
distribution for 3 — 6% central copper-copper col- 
lisions with 35 — 40% gold-gold collisions; the two 



datasets have almost identical N, 



part • 



At least for 



identical nuclei, Npart determines the collision dy- 
namics. It can be used to model other systems, 
such as collisions involving nitrogen. 

The final state composition is well described by 
a thermal model, with production at an equilib- 
rium temperature of = 165 ± 10 MeV; The 
abundance of particle species x, N^ depends on 
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Figure 5. Strangeness suppression factor 7^ as a 
function of Npart- tt, K, p, A, ^ and particles 
and their antiparticles were used to determine 7s 

EQl. 
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its mass rUx [20] : 
Nx: « exp{-m.j,/kTc). 



In contrast to pp collisions, strangeness is fully 
equilibrated; there is no strangeness suppres- 
sion. The system is described by a grand- 
canonical ensemble; a large thermal bath con- 
serves strangeness, and individual strange parti- 
cles can be produced, rather than the pairs re- 
quired in smaller systems. Figure [S] shows how 
the strangess suppression factor 7^ varies with 
Npart- Here, 7s = 1 corresponds to thermal equi- 
librium, Eq. For pp collisions, 7s ~ 0.54. 
Strangeness production is enhanced over pp col- 
lisions even for relatively small systems, such as 
nitrogen- nitrogen collisions. For heavy systems, 
the strangness content almost doubles. 

A more detailed model of the collision includes 
both thermal energy of collective expansion, i.e. a 
hydrodynamically expanding fireball. The 'blast 
wave' model describes this hydrodynamic expan- 
sion with two parameters: system temperature, T 
and a collective expansion velocity, the veloc- 
ity at the outer edge of the fireball f2T . Different 
velocity profiles can be used, with small effects on 
the final results. Because of the collective expan- 
sion, heavier particles have higher (pt)- Fits to 
the pt spectra of different particles find T k, 106 
MeV, and (/?) « 0.55c. This model has been fit 
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to a large body of RHIC data, with considerable 
success PT{ . It described soft-particle production 
in ion collisions, and might be of use for simple 
air-shower simulations. 

Elliptical flow is another aspect of the hydrody- 
namical behavior of the system. In a non-central 
heavy-ion collision, the overlap region is roughly 
almond shaped. Pressure converts this spatial 
anisotropy into a particle density anisotropy. The 
particle flux depends on the azimuthal (per- 
pendicular to the beam) angle with respect to the 
reaction plane (line between the centers of the two 
nuclei). The dependence is 

dN 

— = l + 2i;2Cos(2(j!)) + ... (3) 
dtp 

The elliptical flow, V2 varies with particle species 
and pt- Additional terms for directed flow (e.g. 
Vi) and a small quadrupole moment are ignored 
here. Hydrodynamic flow is a powerful diagnos- 
tic tool, since it probes particle interactions (e.g. 
pressure) very early in the collision. At low mo- 
menta, pt < 2 GeV, the measured particle flow is 
consistent with hydrodyamic models. The initial- 
state spatial anisotropy is completely converted 
into a particle asymmetry. Nuclear matter acts 
like a nearly perfect fluid. 

Figure [6] shows the elliptic flow per constituent 
quark, U2/n {n — 2 for mesons; n = 3 for 
baryons), as a function of pt per constituent 
quark, for a variety of different particles. Ex- 
cept for pions, all of the baryon and meson data 
lie near a single line, with very similar per-quark 
flow [52]. It appears that partons are flow- 
ing, rather than hadrons. Proposed explanations 
for the higher pion V2 include their low mass, 
and/or feed-down from decays of heavier reso- 
nances. Still, the V2 data shows fairly clearly that 
the early interactions involve partons, rather than 
hadrons. 

A final characteristic of the overall event is the 
system size at thermal freezeout. This is deter- 
mined using Hanbury-Brown Twiss (HBT) in- 
terferometry, a measurement of the system size 
that relies on the enhancement in boson (tt, 
etc.) production at small momentum differences, 
\pi — P2I < h/R. By measuring the increased par- 
ticle production at small momentum differences. 
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Figure 6. Elliptic flow per constituent quark, 
V2/n, {n being the number of constituent quarks 
in the hadron) vs. pr/n for various hadrons. 
Except for the tt, the flow of both mesons and 
baryons follows a very similar line, indicating that 
quarks flow, rather than hadrons [22] . 



the source size can be inferred. At freezout (last 
interaction), central gold-gold collisions have a 
Gaussian radius of about 6 fm, about twice the 
size of the original system. 

4.2. Perturbative Probes of the QGP 

Other studies of the QGP use high px probes 
produced in the collisions. High pr particle 
production is well described by pQCD. Except 
for some relatively minor nuclear effects, high 
Pt (usually pr > 2 GeV/c) particle production 
should be the same in pp, pA and AA collisions, 
described by initial state parton distributions, 
pQCD, and universal fragmentation functions. 
Any large differences between systems should be 
due to the interactions between the produced par- 
ticles and the nuclear medium. The time scale 
for fragmentation (whereby partons fragment into 
hadronic jets) is relatively long compared to the 
time a parton remains in the fireball, so medium 
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Figure 7. Raa for tt , 77 and direct photons, to- 
gether with a theoretical calculation of energy loss 
for a plasma with a gluon density of 1100 glu- 
ons/unit of rapidity [M]. Hadrons with a color 
charge are heavily suppressed; photons are not. 



interactions should involve the produced parton, 
rather than the final-state hadrons. Energy loss 
in the medium will manifest itself as a reduction 
in high px particle production as energetic par- 
ticles are shifted to lower pr- This reduction is 
measured by comparing spectra from central AA 
collisions with appropriately scaled pp and/or pe- 
ripheral A A collisions, using the ratio 



In the absence of nuclear effects, Raa = 1. Fig- 
ure [7] shows Raa for tt", t] and direct 7. For direct 
photons Raa ~ 1; photons do not interact with 
the medium. However, for both types of hadrons, 
Raa « 0.2 for 2 < pr < 20 GeV/c, showing a 
large nuclear suppression. In dAu collisions (not 
shown here), Raa ~> 1, as expected due to ini- 
tial state parton scattering. 

In perturbative QCD, the suppression depends 
on the parton density in the fireball and the 
parton-parton cross section. For the standard 
pQCD cross sections, the observed energy loss re- 
quires a parton density of 1100-1200 gluons/unit 
of rapidity 25J. This is far higher than the den- 
sity inferred from the final state multiplicity. 



Raa lias also been measured for heavy quarks. 
Since heavy quarks have a lower velocity (for a 
given momentum) than lighter quarks, less radia- 
tive energy loss is expected. However, Raa for 
heavy quarks is similar to that for light quarks 
|26j . This is difficult to understand in any pQCD 
calculation. 

Raa has also been measured for the J/tp- For 
central AuAu collisions, Raa ~ 0.3 [17] ■ This 
is integrated over all pr, so is not directly com- 
parable to the other Raa measurements. The 
suppression is comparable to what was seen at 
^/Wn = 17.3 GeV at the CERN SPS; this en- 
ergy independence is somewhat surprising. 

Parton energy loss can also be studied with par- 
ticle correlations. A high px particle is selected 
as a 'trigger' particle, and the azimuthal correla- 
tion with a lower px 'associated' particle studied. 
Figure [5] compares the azimuthal correlations for 
dA and for mid-peripheral and central AuAu col- 
lisions, for trigger particles with pT > 8 GeV/c 
and different px associated particles. The dAu 
data has a narrow peak for A(f> ^ (near-side) 
and a slightly broader peak for A<j) « tt (far-side). 
The former is from same-jet correlations, while 
the latter is from correlations involving two back- 
to-back jets; the pp data (not shown) has a simi- 
lar structure. The mid-peripheral AuAu data has 
a similar behavior, although with higher back- 
grounds. In central AuAu collisions, the back-to- 
back peak is much smaller, especialy for softer as- 
sociated particles. The near-side (A(/) « 0) peak 
is largely unchanged. For softer trigger particles, 
the far-side peak largely disappears. 

This data fits a model where the high px parti- 
cles come from parton interactions near the sur- 
face of the fireball. Partons produced deeper in 
the fireball lose most of their energy by dE/dx, 
and so do not produce high px hadrons. High px 
partons produced near the surface, pointing out- 
ward produce full-scale jets, with 'normal' near- 
side correlations. The suppression factor Raa de- 
pends on the surface to volume ratio of the sys- 
tem. Back-to-back jets occur only for a narrow 
range of geometries, where a near-surface inter- 
action produces back-to-back partons nearly tan- 
gential to the surface. 
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Figure 8. Azimuthal correlations between a trig- 
ger particle with pr > 8 GeV/c and an associated 
particle in different pr ranges (in GeV/c) for dAu 
and mid-peripheral and central AuAu collisions 
[281. 



5. Recent Theoretical Developments 

The blast wave model does an excellent job of 
modelling soft particle production at RHIC. How- 
ever, the fit parameters imply very high interac- 
tion cross sections; these cross sections do not 
agree with expectations from pQCD based parti- 
cle interactions. Several related puzzles are also 
of interest: 

1) The elliptic flow is at the hydrodynamic 
limit, with the produced matter behaving like an 
almost perfect fluid. 

2) The Raa measured at high px can only be 
explained in a pQCD framework by an unphysi- 
cally high local parton density. 

3) The Raa for heavy quarks is similar to that 
for lighter particles. 

These puzzles have led to considerable theo- 
retical speculation. Ed Shuryak has proposed 
that the temperature range Tc < T < AT^ is 
a strongly coupled regime for partons [SS]; In 



this non-perturbative phase, there are many very 
weakly bound colored states/resonances, such as 
55 1 9957 etc. These states are very lightly 
bound, so have large radii, leading to rescattering 
cross sections 10 to 100 times those predicted by 
pQCD. Similar behavior is seen for atoms that 
have been tuned (via a magnetic held) to be 
barely bound. Shuryak also pointed out that this 
strongly coupled QGP (sQGP) may be expected 
based on duality arguments with weakly coupled 
string theory. At the same time, lattice calcula- 
tions indicate that meson bound states (notably 
including the J surrvive up to temperatures 
considerably above « 1.5Tc; this could explain 
the lack of additional suppression. 

These strong interactions may explain the large 
elliptic flow and reduced Raa- The strong in- 
teractions could also reduce the QGP lifetime 
(as is seen in some HBT studies). In short, 
this strongly-coupled phase could explain many 
of these puzzling observables. Of course, detailed 
quantitative studies are needed. 

One measure of the cross section is the fluid 
viscosity; elliptic flow depends on the shear vis- 
coity/entropy {ij/s). Flow data shows 77/s < 0.1 
[50] , indicating that the sQGP is a 100 times bet- 
ter fluid than water. This viscosity is much lower 
than that expected for a hadron gas or from a 
perturbative (i.e. hotter) QGP, and appears be- 
yond the reach of perturbative QCD. In fact, it 
approaches the quantum limit, 77/s ~ l/47r, cal- 
culated using duality arguments ^30] . 

6. Lessons for Cosmic-Ray Modelling 

RHIC data offers some guidance for modelling 
cosmic-ray air showers. Nuclear effects are signif- 
icant; A A collisions are not merely superpositions 
of pp or even pA collisions. A blast-wave model 
does a good job of explaining soft particle produc- 
tion; it may be of interest for simple simulations. 
Three aspects of the RHIC data that may be par- 
ticularly significant are: 

1) The reduced forward particle production in 
dA collisions may be due to decreased parton den- 
sities in nuclei at low x. This may reduce number 
of high energy muons, and/or affect the down- 
ward energy fiow in air showers. More quantita- 
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tive saturation models are needed to model higher 
energy interactions. 

2) The strangenss content in nuclear collisions 
is significantly increased over pp collisions even at 
moderate Npart, inpA and AA, and, probably, in 
ttA collisions. At RHIC energies, the increase is 
about 50% for nitogren- nitrogen collisions. Since 

decay faster than tt*, in air showers, they are 
more likely to decay before interacting, so the in- 
creased strangeness should lead to more high en- 
ergy muons. The increased strangeness will also 
reduce the ir'^ fraction in collisions, slowing the 
conversion of hadronic energy into electromag- 
netic energy. This may partially counterbalance 
the previous item. 

3) Lighter systems may be modelled by using 
the Npart dependence of heavy-ion systems. 

7. Conclusions and Future Prospects 

A simple explosive-exansion model can explain 
much of the soft particle production data. Sim- 
ilarly, a surface emission picture can explain a 
lot (but not all) of the high px particle produc- 
tion. However, pQCD calculations do not repro- 
duce the parameters required for these models. 
It may be that a new non-pcrturbative QGP is 
being produced at RHIC energies. 

Over the next few years, the RHIC detectors 
will upgrade their subsystems. The major exper- 
imental goals are high-statustics studies of open 
charm and T production and of leptonic decays of 
vector mesons; the latter may be sensitive to chi- 
ral symmetry breaking. Another interesting idea 
is to search for the tricritical point of QCD; this 
might be found at non-zero baryon density; RHIC 
can search for it point by scanning the machine 
energy. 

I thank the organizers for an enjoyable meeting, 
and my RHIC colleagues for useful suggestions. 
This work was supported by the U.S. DOE under 
contract number DE-AC03076SF00098. 
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